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An algorithm for estimating the optimum combining weights for the Ka-band 
(33.7-GHz) array feed compensation system has been developed and analyzed. The 
input signal is assumed to be broadband radiation of thermal origin, generated by 
a distant radio source. Currently, seven video converters operating in conjunction 
with the real-time correlator are used to obtain these weight estimates. The algo- 
rithm described here requires only simple operations that can be implemented on a 
PC-based combining system, greatly reducing the amount of hardware. Therefore, 
system reliability and portability will be improved . 


I. Introduction 

At the present time, there is considerable interest in 
operating the DSN at increasingly higher carrier frequen- 
cies in order to realize the inherent advantages associ- 
ated with shorter wavelengths, namely, greater antenna 
gains, increased useful bandwidth, and reduced sensitiv- 
ity to plasma effects. Consequently, there is an effort un- 
derway to demonstrate the feasibility of using Ka-band 
(33.7-GHz) carrier frequencies for deep space telemetry. 
However, there are also disadvantages associated with the 
use of higher carrier frequencies, such as greater sensitivity 
to weather effects, increased requirements on pointing ac- 
curacy, and reduced antenna gains due to imperfections in 
the antenna’s reflecting surfaces. Such imperfections be- 
come particularly noticeable on large receiving antennas, 


where gravitational distortions, wind-induced vibrations, 
and collimation problems can seriously degrade antenna 
performance. Some of these losses can be recovered with a 
properly designed compensation system employing an ar- 
ray of receiving horns in the focal plane of the antenna; a 
conceptual design of such a combining system is shown in 
Fig. 1. Complete descriptions and analyses of a real-time 
array feed compensation system designed for deep space 
telemetry can be found in the literature [1—3]. Recently, 
a seven-element array feed compensation system has been 
installed at DSS 13 for the purpose of demonstrating com- 
bining concepts in real time. 

Perhaps the most serious problem encountered during 
the Ka-band array feed compensation effort was the lack of 
reliable coherent sources in the antenna’s far field. Since 


this is not an operational frequency band, no spacecraft 
has yet been built employing Ka-band carrier frequencies 
(Mars Observer did carry a low-power Ka-band beacon, 
but this spacecraft ceased to function before reaching its 
target). However, since only spatial coherence is required 
to demonstrate the array feed combining concept, it is 
possible to carry out weighted combining operations using 
natural radio sources, such as quasars and planets. Since 
these are always in the antenna’s far field, the only addi- 
tional requirement is that the sources remain unresolved. 
This article is devoted to the derivation of the optimum 
combining weights for maximizing he signal-to-noise ratio 
(SNR) of the combined signal and the estimation of the op- 
timum combining weights in real time from natural radio 
signals observed in the presence of additive noise. 

II. The Received Signal 

A functional block diagram of the combining system is 
shown in Fig. 2. The array horns receive a broadband sig- 
nal of thermal origin from a distant point source, each with 
different amplitude and phase introduced by the antenna 
distortions. An independent noise waveform is added to 
each signal in every channel, the result of receiver noise 
plus background radiation received from all directions in 
space. After passing through narrowband filters of center 
frequency w, the received signal can be represented in the 
ifcth channel as 

r k(t) — + n k{t) k — lj 2, ... , A (la) 

where 

s k (t) — y/2S k {a c (t) cos (ut + B k ) 

+ a 3 (t) sin (u)t + #*)} (lb) 

n k (t) = V2{n ck (t)cos(L>t) + n sk (t) sin (u/f)} (lc) 

Note that the signal components in the various chan- 
nels differ from each other only in their amplitude and 
phase, having been generated by the same point source. 
Thus, the random processes s k (t) are completely corre- 
lated. However, the noise processes n k (t) are assumed 
to be uncorrelated in all channels, as these are composed 
of thermal noise generated within the receivers and back- 
ground radiation arriving from all directions in space. 

The received waveforms are downconverted to base- 
band in-phase and quadrature signals rj k (t) and rq k (t) 


by premultiplying with local oscillator signals of the form 
\/2cos (u>t) and \/2sin (w<), and low-pass filtering: 


rn(t) = «/fc(0 + n Ik (t) (2a) 

r Qfc(0 = s Qk(t) + «QJb(0 (2b) 

s fk {t) = S k {a c (t) cos (0 k )-a s (t) sin(0*)} (2c) 

SQ k {t) = S k {a c (i) sin (0*) + a,(<) cos(0*)} (2d) 

nj k (t) - n ck (t) - n sk (t) (2e) 

n Qk (t) = n ck (t) + n sk (t) (2f) 


After sampling the baseband waveforms, the resulting 
in-phase and quadrature samples may be treated as com- 
plex samples rjt(z), defined as 

r k (i) = s k (i) + n k (i) (3a) 

where 

h(i) = s Ik (i) + js Qk (i) (3b) 

n k {i) = n Ik (i)+ jn Qk (i) - (3c) 

In other words, we shall use complex arithmetic to operate 
on these samples. 

Defining the complex coefficient S k = S k ^ 6k , Eq. (3b) 
can also be written as 

s*(i) = § k {a c {i) + ja s (i)} = S k a(i) (4) 

which shows that the complex channel scaling factors sep- 
arate from the temporal variations. 

The real and imaginary parts of the complex noise se- 
quence are independent random sequences, each with vari- 
ance <J k . The components of the signal sequence, a c (z) 
and a 3 (z), are also independent random sequences, being 
of thermal origin. However, since we are interested in ex- 
tracting the complex magnitude S kt it is reasonable to let 
each signal component have variance 1/2, so as not to in- 
troduce additional scaling. In the following, the complex 
signal coefficients and the noise variances will be assumed 
to be constants independent of time. 
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III. Combining Weights to Maximize SNR 

The goal of the combining operation is to maximize 
the SNR of the combined sequence. The approach is to 
multiply each sequence by a complex combining weight 
with the property that the sum of the weighted sequences 
achieves the greatest SNR. 

Let the Jcth combining weight be denoted by w k - Mul- 
tiplying each received sequence by the corresponding com- 
plex weight and summing yields the combined sequence 
*(*): 


E (VWO (?fe-) 

E2KI 2 <t| 


K 





A po 


(9) 


which shows that for any choice of combining weights, the 
achievable SNR is bounded from above by the sum of the 
channel SNRs. Except for an arbitrary complex factor, 
equality is achieved when we let y/2w\ ter* = Sl/V2a ki 
whereby the optimum combining weights are determined 
in terms of the signal and noise parameters as 


K 

z(0 = 53 a (0 

k - 1 


(5) 



( 10 ) 


With s c (i) = Ef=i^(*) **(»') denoting the combined 
signal, the SNR of the combined sequence is defined as 


_ E\s c (i)\ 2 

P* var{z(z)} 


k—i 


E 2|«*| 3 *2 

k = 1 


( 6 ) 


This is the quantity we wish to maximize by judicious 
selection of the combining weights. 

The optimum combining weights are obtained by means 
of the Schwarz inequality. Writing the combined sequence 
as 


The combined SNR is maximized when each sequence 
is multiplied by the conjugate of the signal scaling factor 
and divided by the variance of the additive noise in that 
channel. Therefore, these quantities have to be estimated 
in real time to determine the correct combining weights. 


IV. Parameter Estimates 

The estimator described here is based on the observa- 
tion that the temporal variation of the signal components 
is identical in every channel. This implies that the ex- 
pected value of the product of a received sequence with the 
complex conjugate of a sequence from any other channel 
is equal to the product of the complex coefficients. That 
is, 


t = i 



(7) 


E [r;(i)»Vn(*)] = E [(&a(t) + «/(*)) (S m 5(0 + «m(*))] 

= sis m ( 11 ) 


and applying the Schwarz inequality, yields 


y (yzwkVk) 

k - 1 



K 


y2\w k f 


Kk = l 




J 


The last equality follows from the definition E\a(i )\ 2 = 
1, and the assumption that noise sequences are not cor- 
related with other noise sequences or with the signal se- 
quence. If the received sequences are ergodic, then ensem- 
ble and time averages are identical, suggesting the follow- 
(8) ing estimates for the coefficient products: 

Si Sm 1 / *\ 


Dividing both sides by the first term on the right-hand 
side of Eq. (8), we obtain 


t,m= 1,2, ... ,K; l±m (12) 
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As the number of terms grows without bound, the es- 
timation error approaches zero. With the total number 
of channels equal to K , index the channels according to 
the SNR, so that channel number 1 contains the greatest 
SNR, channel number 2 the second-greatest, and so on (in 
case of equalities, an arbitrary choice can be made). Thus, 
we may view channel 1 as the “reference channel,” as it 
provides the highest fidelity signal. Consider the estimates 
of the products S{S m • 

— 1 1 

= m = (13) 

i= 1 

If were known, then we could estimate 5 m , m = 
2 , 3, . . . , K, by means of the formula 


and 


E[Pm 


= 2a 


2 

1 


(16b) 


it follows that 



Ptl — Pnl\ Ptl > Pn 1 
not defined; P ni > P tl 


(17) 


is a reasonable estimate of |5i| 2 when the average signal 
power is large compared to the random variations. When 
the signal power is not sufficiently great, it is possible for 
the estimate to become negative, which is meaningless for 
a power estimate and, therefore, should not be used. 


S{Srn 

Si 


(14) 


The estimate of 


Si 


follows from Eq. (17) as 


Even if only \S\\ were known (that is, if no phase in- 
formation could be obtained), we could still determine 
where 0 1 is the argument of Si. Since in our 
application the combining weights may be multiplied by 
a complex constant without affecting the combined SNR, 
the phase of the reference channel need not be estimated. 
Thus, we turn our attention to the estimation of the signal 
magnitude in the central reference channel. 

An estimate of |5i[ can be obtained in the following 
manner, using well-established experimental techniques. 
With the antenna pointing “on-source,” we obtain an es- 
timate of the total power Pt\, by averaging L independent 
sample powers: 


c i J; . 1 2 

Pn — 7 / , Sia(i) + ni(i) 


»=i 


(15a) 



(18) 


Using this estimate, the signal coefficients for the outer 
channels may be obtained from Eqs. (12) and (14), except 
for a common complex coefficient, as follows: 


-W i 


Q* Q 

El 


(19) 


The complex coefficient e~ i 01 has no effect on the com- 
biner performance, since the combining weights will in- 
clude this factor, effectively setting the average phase of 
the combined signal to zero. 


Next, we point the antenna “off-source” and obtain an 
estimate of the noise power using subsequent samples: 


2 L 


Pnl = J 


i=L+ 1 


(15b) 


Since 


V. Statistics of the Weight Estimates 

The combining weights that maximize the SNR of the 
combined samples are defined in Eq. (10), repeated here 
for convenience: 


Wk 



(20) 


E 



Each complex weight may be multiplied by a complex 

(16a) 

constant without affecting combiner performance. Thus, 


&9 



for our purposes, it is sufficient to estimate the “rotated” 
weights 


wic e ~ j91 


S' k e-j» _ S^_ 

2 2| SiWl 


( 21 ) 


instead, without explicitly determining the phase of the 
central channel, 0\. 

The denominator of each weight is simply an estimate 
of the noise variance in that channel. The noise variance 
in the kth channel, 2 <j\, is estimated by means of simul- 
taneous off-source power measurements in all K channels 
(both real and imaginary noise components have variance 
(t\ ) hence the factor of 2 in the variance of the complex 
noise). If M samples are used, the estimate of the variance 
in the Arth channel is of the form 


product Si§l (whose expectation is proportional to the 
correlation between the first and kth channels), this ap- 
proach also requires an estimate of the magnitude of the 
complex scale factor in the central channel. This quan- 
tity can be updated with each new weight estimate, us- 
ing the last noise power measurement, or a separate on- 
source measurement can be made periodically for the cen- 
tral channel as well. Both cases will be considered. First, 
assume estimates are made with each update, using the 
same number of samples as for the weights (L). The esti- 
mate of the signal power in the central channel then be- 
comes 


\Si 


L + M 


= X^|^ l5 (') + fi 1 (*)| “jjj E l"i(‘)| 2 (25) 

i = 1 t = L + l 


with expectation 


_ i m 

^ = 2M E 

( 22 ) 

, r m m 'j 

= + E n /(*)j 

Clearly, the expected value of this estimate is the de- 
sired quantity, 

= * l ( 23 ) 


E 




(26) 


and variance 





(27) 


Again, we observe from Eq. (26) that an unbiased es- 
timate is obtained. If a separate power measurement is 
carried out for the central channel using N independent 
samples, then the signal power estimate becomes 


hence the estimator is unbiased. Assuming the underlying 
processes are Gaussian, the estimate is a scaled central 
chi-square random variable with 2M degrees of freedom 
and variance 



i J 2 1 ^ 

i£|s,«(;) + a,(0 -jj Y, l»i(OI ! (28) 


V “W) = 4^2 < 4 "^) = M (24) 

In typical applications, the component of the noise vari- 
ance due to the receiver remains constant, while the com- 
ponent due to the background changes slowly with ele- 
vation and azimuth, but may be considered constant for 
many weight estimates. The frequency with which off- 
source measurements must be performed depends on the 
specific details of the experiment. 

The numerator of each rotated weight estimate can be 
obtained using Eq. (19). In addition to an estimate of the 


This estimate is also unbiased; hence, Eq. (26) still 
holds. However, the variance of the estimate decreases 
if N > L, as shown by comparing the following expression 
with Eq. (27): 


var^ 


(KK( 


2 ** + 



(29) 


The estimate of the signal power can, therefore, be 
made arbitrarily good by making both N and M suffi- 
ciently large, provided the natural time scales of the rele- 
vant process variations are not exceeded. 
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Next consider the statistics of the coefficient product 
estimates, SjS* k . As shown in Eq. (11), the expectation 
of the product of the received samples between any two 
channels is the desired coefficient product, implying the 
unbiased estimator structure defined in Eq. (13). 

The variance of this estimate can be obtained by first 
deriving the second moment of the coefficient product, 
subtracting the square of the expectation, and dividing the 
result by the total number of independent sample products 


ElSjSH 2 = E{ (s,S;|5| 2 + Sjahl + S^nj 

= E {|Sj| 2 |S*| 2 |a| 4 + ISjflaflntl 2 + | 2 |a| 2 | 2 + |r* i | 2 |«*| 2 } 

-f (products whose expectation goes to zero) 

= & | 2 | 5 * I W + I §J ?W l **] 2 + \Sk I 2 H 2 N 2 + \nj i 2 |«*| 2 ( 31 ) 


averaged to obtain the estimate. If L sample products are 
averaged, this can be expressed as 

var = 1 j# \SjS* k \ - ^E (?,£*) | ) (30) 

Consider the second moment first. Writing the received 
samples as in Eqs. (3a) and (4), the second moment of the 
sample products between two distinct channels j and k 
becomes 


hjnt ) (s*S*|a| 2 + s;a*h k + S k ah * + njn*) } 


where the overbar denotes expectation. Letting \hj\ 2 = 
2 ( 7 ? , jap = 1, and with \a\ 4 = 2, the second moment 
becomes 

£|5j^| 2 = 2|5j| 2 |5jfc| 2 + 2)S,| 2 <r 2 + 2|S*|V? + a 2 k a 2 

(32) 

Subtracting the square of the expected value and divid- 
ing by L yields an expression for the variance: 

var (SjSlj = \Sj\ 2 \S k \ 2 

4 2 (\Sj\ 2 <?l 4* |5jb| 2 <r?^ + 4(7 ?(T 2 (33) 

If the estimation errors are sufficiently small, each esti- 
mate may be written as its true value plus a small random 
deviation. Assuming this to be the case, we have 


SiS k = S* 1 S k ( l + a) 

(34a) 

|li|2 = |5i| 2 (l + /?) 

(34b) 

=o- 2 (! + t) 

(34c) 

Using Eqs. (24), (27), and (33), it follows that each 
error term is a zero-mean random variable with variance 

, ^ 2 
var (a) = — 

L 


x / 1+ +2 ( + a " ) 
l |5i| a |5 t | 3 \\S k \ 2 |5i| 2 J 

) (35a) 


var (/?) =]i^ { h ( 2<r? + i51 12 ) + ^ } ; 

Q = Lot N (35b) 
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var ^ = M 


(35c) 



(41) 


The estimate of the rotated weight can now be written 
as 


and 


SiSj (l + a) 

2*2 & I y/fjL+mi + y) 


var (() = E 


a 


p p ^ p l 2 

— — 7 — a— — 07 + 7 — 


~ (1 + Q ) 

‘(l + f) (1 + 7) 


(36) 


— B 2 B 2 B 2 

= a 2 + — -f 7 2 + 7 2 — + » 2 7 2 + “J’T 2 

+ (products whose expectation is zero) (42) 


where Ck AS\SI/2 <t\\S\ |. Expanding the denominator, 
multiplying through, and keeping terms up to second order 
yields 


= Ck 



1 & 

l--+a-7 


+ a- 7 (l + a- 7 ) 


(37) 


The expected value of this expression is 


E(w k e -i«>) 



Since i?(a) = £(/?) = £(7) = 0, and a, /?, and 7 
are statistically independent, it follows that all terms in 
Eq. (42) containing a , /?, or 7 to the first power go to 
zero. Thus, the normalized variance of the rotated weight 
may be expressed in symmetric form as 


\Ck\ 2 var ^w>jfce— ) 


var (a) + var 


(S) 


+ var (a) var ^) 
+ var (ct) var (7) 


4- var (7) 


+ var (7) var 



(43) 


= C k 


1 + var 



+ var (7) 


(38) 


The statistics of these rotated weight estimates for spe- 
cific parameter values are examined in the next section. 


which shows that the estimate is biased, but becomes 
asymptotically unbiased as the number of samples grows 
without bound. Writing the rotated weight estimate as 

= £*(1 +0 ( 39 ) 

where £ = a — ^ — 7 — 07 — «f-h7f + ^- + 7 2 , it follows 
that 

var = |Cfc| 2 var (£) (40) 

Here £ is a random variable with second-order statistics 


VI. Numerical Results 

The expressions for the normalized bias and variance 
of the rotated weights developed in Eqs. (41) and (42) 
are functions of the parameters A/, N, and L } as well as 
of the signal powers and noise variances in the combiner 
channels. The noise variance is largely determined by the 
physical characteristics of the front-end amplifiers, and in- 
cludes contributions from background radiation as well. 
The signal power in each channel depends on the strength 
of the source, the antenna aperture, and the amount of 
distortion suffered by the main reflector. However, the re- 
maining parameters can be selected to achieve a desired 
level of estimator performance. 
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In a typical combining experiment, the antenna is first 
pointed off source, and the noise power in each complex 
channel is estimated according to Eq. (22). Denote the 
number of samples used to make this estimate M. Next, 
the antenna is pointed on source, and the signal in the cen- 
tral channel is measured by means of Eq. (25) or Eq. (28), 
using either L or N samples, where L is the number of sam- 
ples per update. Of course, the signal levels in the other 
channels could also be determined, but that information is 
not used by this algorithm; instead, the rotated complex 
weights are estimated according to Eq. (19). Theoretically, 
some improvement could be obtained if independent signal 
power measurements were used to improve the magnitude 
estimates, but that issue will not be addressed here. 

After the coefficients are determined, combining weight 
updates are obtained every L samples for use in the real- 
time combining operation. The inequality (M, N) » L is 
usually valid, since both signal and noise power levels tend 
to remain constant over a great many coefficient updates. 
However, we shall also consider the case N = L, corre- 
sponding to a situation where the signal power in the cen- 
tral channel changes fast enough to warrant its measure- 
ment with each update. This could occur if wind-induced 
dynamics or time- varying pointing errors were present. 

Representative values were chosen for the signal power 
levels, normalized to the noise variance, which was as- 
signed a value of 1 power unit. With a channel noise tem- 
perature of roughly 100 K, a strong Ka-band radio source 
(such as Venus or Jupiter) may produce a 10-K rise in the 
central channel when observed with the 34-m antenna at 
DSS 13. At low and high elevations, gravitationally in- 
duced distortions of the main reflector generally deflect as 
much as 10 percent of the total signal power to some of 
the outer channels. Thus, the values — l,|Si| 2 = 0.1 
and |5jk| 2 = 0.01, k ^ 1 will be assumed for the numerical 
examples that follow, meaning that the signal power in the 
central (reference) channel is one-tenth of the noise power, 
while typical signal powers in the outer channels are 10 
times smaller than those in the central channel. 

The normalized bias defined in Eq. (41) is shown in 
Fig. 3 as a function of A, with M as a parameter. For any 
AT, the minimum achievable bias is a monotone decreasing 
function of M that approaches an asymptotic limit from 
above. Minimum values of M and N are clearly specified 
for any desired bias level: for example, if we wish to main- 
tain a normalized bias of 0.001, then M must be greater 
than about 10 5 and N must be at least 3 x 10 5 ; at a sam- 
pling rate of 2 x 10 5 samples/sec, this would take a mere 


2 x 1.5 = 3 sec, which is short compared to typical time 
scales encountered in practice. If we let M exceed 10 6 , 
then the above value of N will suffice. Suppose we select 
this value of N in order to meet the bias requirements, and 
examine the normalized variance defined in Eq. (43). 


The normalized variance of the rotated weight esti- 
mates is shown in Fig. 4 as a function of L. It is clear 
that a given value of N specifies a limit for the smallest 
attainable variance and, hence, limits the performance of 
the weight estimator. For the specified value of V, this 
limit is about 4 x 10“ 4 , even as M and L grow without 
bound. The dependence on N disappears if the signal 
power is recomputed with each update, using the previous 
L samples, corresponding to the case N = L. Now any de- 
sired level of performance can be achieved provided that 
both M and L are sufficiently great, as shown in Fig. 5. 
Given M and L, the corresponding normalized bias can be 
determined from Fig. 3 by substituting L for N . Note that 
for very small values of L, the variance increases rapidly 
when N = L> due to the greater error in the signal power 
estimate as a result of insufficient observations. 


VII. Summary and Conclusions 

A digital combining-weight estimation algorithm for use 
with broadband sources has been described and analyzed. 
Although the algorithm provides a biased estimate of the 
combining weights, the bias can be reduced to any de- 
sired level by observing enough samples. The normalized 
variance of the weight estimates can be similarly reduced; 
however, care must be taken to obtain accurate signal and 
noise power estimates, as these quantities ultimately limit 
the performance of the estimator at all update rates. With 
this algorithm, only three primitive estimates are needed 
to obtain each combining weight estimate: the noise vari- 
ance in the kth channel, the signal power in the central 
channel, and the complex correlation coefficient between 
the central and fcth channels. The first two are simple 
power measurements requiring magnitude squaring and 
accumulation, while the third is a complex multiply-and- 
accumulate operation; both operations can be carried out 
with current PC-based digital hardware, provided the sam- 
pling rates do not exceed roughly 250,000 samples/sec. 
Consequently, this approach is useful for reducing estima- 
tor complexity while greatly increasing the portability of 
the entire array feed compensation system, enabling the 
demonstration of combining gain on a variety of DSN an- 
tennas. 
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Fig. 1 . Real-time antenna-compensation system conceptual design. 
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Fig. 2. Signal-combining system (baseband model). 
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